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Abstract 

A strict positivity of the ground-state energy is a necessary and sufficient condition 
for spontaneous supersymmetry breaking. This ground-state energy may be directly 
determined from the expectation value of the Hamiltonian in the functional integral, 
defined with an antiperiodic temporal boundary condition for all fermionic variables. 
We propose to use this fact to observe the dynamical spontaneous supersymmetry 
breaking in Euclidean lattice simulations. If a lattice formulation possesses a manifestly 
preserved fermionic symmetry, there exists a natural choice of a Hamiltonian operator 
that is consistent with a topological nature of the Witten index. We numerically 
confirm the validity of our idea in models of supersymmetric quantum mechanics. 
We further examine the possibility of dynamical supersymmetry breaking in the two- 
dimensional A/" = (2,2) super Yang-Mills theory with the gauge group SU{2), for 
which the Witten index is unknown. Although statistical errors are still large, we 
do not observe positive ground-state energy, at least within one standard deviation. 
This prompts us to draw a different conclusion from a recent conjectural claim that 
supersymmetry is dynamically broken in this system. 
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§1. Introduction 

The possibility of the spontaneous breaking of supersymmetry (assuming that it is not 
broken at the tree level) is a highly dynamical issue and its precise study requires a nonper- 
turbative framework. Generally, the Witten index^^ provides an important clue. One can 
infer that dynamical supersymmetry breaking does not occur in a wide class of supersym- 
metric models where the Witten index can be computed to be nonzero. However, the Witten 
index is not a panacea. There are still many interesting models for which it is very difficult to 
determine the Witten index and, in some cases, the index itself would be ill-defined because 
of a gapless continuous spectrum.^)' ^^ 

On the other hand, it is well known that a strict positivity of the ground-state (or vacuum) 
energy is a necessary and sufficient condition for spontaneous supersymmetry breaking.^) In 
principle, therefore, one can judge whether supersymmetry breaking occurs by computing 
the ground-state energy. 

In this paper, in light of recent developments on the lattice formulation of supersymmet- 
ric theories, ^)"^) we propose to observe dynamical supersymmetry breaking with Euclidean 
lattice formulations employing the above idea. This work was originally motivated by a re- 
cent paper by Hori and Tong^^ in which they conjectured dynamical supersymmetry breaking 
in the two-dimensional M = (2,2) super Yang-Mills theory with the gauge group SU{Nc). 
Lattice formulations of this two-dimensional theory are the simplest among recent lattice 
formulations of extended supersymmetric gauge theories. ^°-'"^^-' Therefore, it is highly natu- 
ral, if it is possible, to examine supersymmetry breaking in the two-dimensional A/" = (2, 2) 
super Yang-Mills theory with lattice formulation. *) This is what we do in this study. To 
our knowledge, this is the first instance in which the dynamical supersymmetry breaking 
in gauge field theory is investigated numerically, although there exists a closely related 
and thought-provoking observation in Ref. 19). (There are a number of numerical works 
related to this issue in one-dimensional supersymmetric models^^-*"^^-* and two-dimensional 
Wess-Zumino models. ^^•'"^^•') Although statistical errors in our Monte Carlo study using a 
formulation in Ref. 13) are still large, we do not observe positive ground-state energy, at 
least within one standard deviation. This observation prompts us to draw a different con- 
clusion from the conjectural claim in Ref. 9). This is the content of §4. In §2, we present 
our basic idea concerning the determination of the ground-state energy in the Euclidean 
functional integral formalism. Then, in §3, we illustrate how our method works by applying 
it to supersymmetric quantum mechanical models.^-* Section 5 is devoted to the discussion. 



*' While preparing this paper, we discovered a preprint^^-* in which this problem is addressed on the 
basis of the lattice formulation in Ref. 11). 



In what follows, the boundary condition of fermionic variables for the temporal direction, 
whether it is periodic (PBC) or antiperiodic (aPBC), is crucial. For all bosonic variables and 
for all variables with respect to the spatial directions, we will assume the periodic boundary 
conditions. Unless noted otherwise, the term "boundary condition" will always refer to the 
boundary condition of fermionic variables for the temporal direction. 

§2. Basic idea 

What we want to determine is the ground-state energy Eq of supersymmetric theories. 
If Eq > 0, supersymmetry is spontaneously broken and it is not if Eq = 0.^-* With the 
Euclidean functional integral formalism, one could determine the ground-state energy from 
the expectation value of Hamiltonian H, 

{H}pBC-—r -. Z^, l^-ij 

JpBc ^^ ^ 

where S is the Euclidean action and dyU symbolically denotes a measure for the functional 
integration. We assumed the periodic boundary condition (PBC) for fermionic variables 
because this boundary condition is consistent with supersymmetry. One would then be able 
to obtain the ground-state energy Eq by taking the large imaginary-time limit /9 — > oo, 
where /3 is the temporal size of the system,*^ because in this limit, only the contribution of 
the ground-state (s) survives in Eq. (2-1). 

However, this naive idea is wrong. First, we must note that the functional integral in the 
denominator of Eq. (2-1) is proportional to the Witten index,^^)'^^)**) 

ZpBc = ATpBc / c^/xe-^ = Tr(-l)^e-'3^ = Tr(-l)^, (2-2) 

where F is the fermion number operator and A/pBc is a proportionality constant that depends 
on the choice of the integration measure dfi. The constant A/pBc may depend on ultraviolet 
and infrared cutoffs (the number of lattice points for lattice regularization) and possibly on 
the boundary condition. Second, the numerator of Eq. (2-1) is proportional to the derivative 
of the Witten index with respect to P, which is always zero: 

AApBc / dfi He-^ = Tr(-l)^iJe-^^ = -^ Tr(-l)^e-^^ = 0. (2-3) 



*' f3 is not to be confused with the conventional gauge coupHng constant in lattice gauge theory. 
**^ Our discussion in this section is based on the assumption that the expressions appearing in Eq. (2-2) 
are meaningful. For this, we may assume that the spectrum of H is discrete so that the Witten index is 
unambiguously defined. Our basic formula (2-7) for the ground-state energy itself, however, might also be 
applicable to systems in which this assumption fails. 



This independence of the Witten index from a parameter of the theory, (3, is a consequence 
of the supersymmetry algebra.^-* Thus we have 

(^)PBC = ;^r^. (2-4) 

We finally recall that the Witten index vanishes when supersymmetry is spontaneously bro- 
ken. Therefore, we see that {H)pbc is indefinite when supersymmetry is broken. Otherwise, 
it is zero or indefinite depending on whether or not the Witten index is nonzero. Note also 
that a similar remark is valid for the expectation value of generic operators when the peri- 
odic boundary condition is imposed; when the Witten index vanishes, ^pbc cannot be used 
as a normalization factor for expectation values. With the periodic boundary condition, 
therefore, the expectation values normalized by the partition function can be ill-defined, and 
in such a case, we must consider "denominator-free" expectation values, such as Eq. (2-3). 
In any case, the expectation value {H)pbc does not provide useful direct information on the 
ground-state energy or on supersymmetry breaking. 

In Eq. (2-3), what prevents us from obtaining the ground-state energy is the factor 
(—1)^, which is the heart of the Witten index. If this factor can be removed, the above idea 
of using the expectation value of the Hamiltonian would be valid. As is well known (see, 
for example, Refs. 39), 40)), such a removal can easily be achieved. What we must do is 
simply to change the boundary condition of fermionic variables for the temporal direction 
from periodic to antiperiodic (aPBC). This defines the thermal partition function with the 
inverse temperature f3, instead of the Witten index, 

Z,pBc = A4pbc / dfi e-' = Tr e"'^^, (2-5) 

JaPBC 

which should be positive definite, and, as its derivative with respect to /9, 

Kpbc f dfiHe-^ = TTHe-'^". (2-6) 

iaPBC 

Therefore, taking the long-time limit (or the low-temperature limit) of the ratio of these two 
quantities, we have 

lim f^ aPBC = lim if^c^ ^ j^^ ^ ^ 2-7 

and ground-state energy E^) is obtained. This is our basic formula. 

Before taking the /9 — ^ oo limit, Eq. (2-7) is merely the expectation value in the thermal 
equilibrium with finite temperature 1//9, and supersymmetry is explicitly broken by the 
temperature. In this aspect, it is interesting to note an analogy to a conventional way of 



detecting the spontaneous breaking of an ordinary symmetry, for example, the Z2 symmetry 
of the Ising spin. In this case, one breaks the symmetry by applying an external magnetic field 
that is conjugate to the order parameter of symmetry breaking, that is the magnetization. 
One then observes (in the thermodynamic limit) how a trace of the breaking remains after 
the applied field is turned off. 

For supersymmetry, the order parameter is a positivity of the ground-state energy and the 
conjugate variable to the energy is the temperature. In Eq. (2-7), we break supersymmetry by 
placing the system in thermal equilibrium. We then observe how the effect of the temperature 
remains in the zero-temperature (or the large imaginary-time) limit /9 — * 00, for which one 
naively would expect that the effect simply disappears. If the effect remains, we judge that 
spontaneous supersymmetry breaking occurs. Recall that the expectation value (-ff)pBC 
with the periodic boundary condition is always zero or indefinite. Thus, if a well-defined 
-Eq > is found through Eq. (2-7), it is the effect of the boundary condition surviving even 
in the long-time (or the zero-temperature) limit. *^ Physically, this survival of the effect of 
the boundary condition can be understood in terms of the appearance of a massless (or 
zero energy for quantum mechanics) Nambu-Goldstone fermion associated with spontaneous 
supersymmetry breaking. 

Our basic formula (2-7) is very simple. However, to embody it in Euclidean lattice 
formulation, there are still several issues to be clarified. First of all, the above argument 
assumes that regularization to define the functional integral does not break supersymmetry. 
Otherwise, one would not be able to distinguish spontaneous supersymmetry breaking from 
a possible explicit breaking due to regularization. As is well recognized, generally, a regu- 
larization based on a spacetime lattice is irreconcilable with supersymmetry. For theories 
with the extended supersymmetry, it is nevertheless sometimes possible to set up a lattice 
regularization that preserves the invariance under some supersymmetry transformations. ^-'"^•* 
Then, if the spacetime dimension is low enough, one may expect that the invariance under 
a full set of supersymmetry transformations is restored in the continuum limit. In what 
follows, we assume this sort of lattice regularization. In particular, for the study of the 
two-dimensional M = (2, 2) super Yang-Mills theory, we adopt the formulation in Ref. 13) in 
which a fermionic symmetry Q that is a part of the supersymmetry is manifestly preserved. 
We will briefly review this formulation in §4. 

Secondly, closely related to the above point, we must properly choose a possible additive 
constant in the Hamiltonian H. In other words, we must correctly choose the origin of the 
energy. This point is, of course, crucial for judging spontaneous supersymmetry breaking 



*' Note, however, that spontaneous supersymmetry breaking differs from spontaneous breaking of ordi- 
nary symmetries in that it can occur even in a system with finite volume.*' 



from the positivity of Eq. Note also that, when the Witten index is nonzero, supersymmetric 
invariant state(s) must have a precisely zero energy eigenvalue Eq = for relation (2-3) to 
hold. That is, Eq. (2-3) is not invariant under an arbitrary shift of the origin of the energy 
H ^ H + c, when the Witten index is nonzero. 

Of course, a natural prescription for defining the Hamiltonian is to use the supersymme- 
try algebra. One may first define supercharge operators Q and Q^ (with some regulariza- 
tion) and define a (regularized) Hamiltonian operator H by the anti-commutation relation 
H = {Q, Q^}/2 without any additive constant. This is precisely the idea behind the Hamil- 
tonian formulation of supersymmetric theories^^-'"^^-''^^-''^^'''^^-''^^-''^^-*"^^-* with which a possible 
additive constant in the Hamiltonian H is automatically fixed. *•* 

For the following reason, however, this issue of a "correct" Hamiltonian is somewhat 
delicate in the functional integral formulation based on the Lagrangian. 

Suppose that the (for simplicity, off-shell) supersymmetry algebra is realized by the trans- 
formation law for variables appearing in the continuum Lagrangian. This implies that there 
exists a fermionic transformation Q such that {Q, Q} = 2ido, where Q is the conjugate 
fermionic transformation of Q and do is the time derivative. One would then expect, from 
this algebra, that the relation iQQ = 2H holds,**-* where Q is the Noether charge***^ as- 
sociated with Q and H is the Hamiltonian obtained from the Lagrangian by the Legendre 
transformation. If this relation holds, this H could be used in the functional integral as a 
Hamiltonian operator that is consistent with the supersymmetry algebra. 

In reality, however, the relation holds only up to equations of motion. Generally, one 
ends up with 

-QQ = H + (terms being proportional to equations of motion). (2-8) 

That is, a Hamiltonian suggested from the algebra can differ from the original one obtained 
through the Legendre transformation from the Lagrangian. This occurs very commonly, 
and we will encounter such a situation even in the simplest supersymmetric system in the 
next section. The additional terms, which would be negligible in the classical level, cannot 
be neglected in general within the functional integral because those terms may give rise to 



*' On the other hand, from the viewpoint of the feasibihty of numerical simulations that preserve the 
gauge symmetry, the Euclidean lattice formulation appears advantageous. 

**) We take normalization of the Noether charge such that the Poisson bracket {Q, -jp generates the Q 
transformation. After the quantization, the relation would be read as H = {Q, Q}/2 which is consistent 
with the positivity of the Hamiltonian. 

***) iiY field theories, when supersymmetry is spontaneously broken, the Noether charge (supercharge) 
itself would be ill-defined owing to a massless singularity associated with the Nambu-Goldstonc fcrmion. In 
the field theory case discussed in §4, we use the Noether current instead. 



contact terms at a coincident point, i.e., ultraviolet-divergent constants. We thus have two 
(among possibly many) options for a Hamiltonian operator in quantum theory: one is the 
original Hamiltonian obtained from the Lagrangian and the other is iQQ/2. How can we be 
sure that we are using a Hamiltonian with a correctly chosen additive constant before we 
measure the ground-state energy? Clearly, we need some guiding principle. 

We have no general answer to the above question. See also §5. However, if the lattice 
formulation one adopts possesses at least one exactly preserved fermionic symmetry, for 
example, the above Q, there exists a natural prescription for a choice of the Hamiltonian. It 
is the left-hand side of Eq. (2-8), H = iQQ/2. This choice of the Hamiltonian in quantum 
theory corresponds to "renormalizing" additional terms in the right-hand side of Eq. (2-8) 
into H. This definition is natural because the structure H = iQQ/2 is suggested from the 
supersymmetry algebra. Moreover, this choice has the correct origin of the energy in the 
sense that it is consistent with the topological property of the Witten index, Eq. (2-3). That 
is, 

AfpBC [ dfiHe'"' = XPBC [ dfi'-QQe-' = XpBc [ dfiQ (^Qe''] = 0, (2-9) 

JPBC JPBC ^ JPBC V^ / 

where we have used the Q-invariance of the action and of the integration measure. *) As 
already noted, when the Witten index is nonzero, this property fixes the origin of the energy 
uniquely. For these reasons, we consider that the definition H = iQQ/2 is natural. Of course, 
for cases of interest, we do not know a priori whether the Witten index is nonzero or not, 
and if it is zero, the above argument for the structure H = iQQ/2 based on relation (2-3) is 
groundless (a shift of the origin H -^ H+c does not influence Eq. (2-3) if A/pBc /pbc ^^^ ^~^ ~ 
0). Nevertheless, we adopt this Q-exactness of the Hamiltonian as a working hypothesis in 
what follows because the definition of a Hamiltonian operator should be independent of 
whether or not the supersymmetry is spontaneously broken. 

§3. Supersymmetric quantum mechanics 

In this section, we examine our method by applying it to a Euclidean lattice formulation 
of the supersymmetric quantum mechanics.^-* We find that this example provides a good 
illustration of our method. 

The Lagrangian of the supersymmetric quantum mechanics is given by 

L = ^(50)^ - ^{Wr + ii{td - W")i, + \f\ (3-1) 



*' Strictly speaking, to show this relation, we must assume that the integral /pg(^ d^Qe is finite. 



where all variables are functions of the "time coordinate" x and d is the derivative with 
respect to x, (9 = d/{dx). and F are bosonic variables and ip and ip are fermionic. The 
superpotential W = W{(f)) is a function of and the prime denotes the derivative with 
respect to 0. With the periodic boundary condition for all variables, the action S = j dx L 
is invariant under the following A/" = 2 supersymmetry transformations: 



and 






w\ 



QF = -idip + W'lp, 



(3-2) 
(3-3) 



^, 



QF = id^ + W"i). 



Q^ = -F + id(p + W, 
One can confirm that the transformations form the supersymmetry algebra 



^t2 



Q' = Q =0, {Q,Q} = 2td 



(34) 
(3-5) 



(3-6) 



off-shell, i.e., without using any equations of motion. In this system, it is well known"^) 
that supersymmetry is spontaneously broken if and only if the number of zeros of the func- 
tion W'{(f)) is even, or equivalently, W^(— oo) and W{+oo) have opposite signs. (We assumed 
that |Vr(±oo)| = +00.) 

A crucial fact for us is that the classical action can be expressed as the Q-exact form*^ 



S= f dxL = Q f dx-ip{F-id(l) + W'). 



(3-7) 



Then the invariance of S under Q and Q can be easily seen by using the supersymmetry 
algebra (3-6).**) 

The Hamiltonian corresponding to the Lagrangian (3-1) is given by 



and, as we noted in Eq. (2-8), we have 

^QQ = H+ -F{F - id<P - W) + \i^{id - W")ij, 

Zi Zi Zi 

where Q, is the Noether charge associated with the Q invariance: 



(3- 



(3-9) 



(3-10) 



*^ For this, we must note that / dx iW'dcj) — J dx idW — 0. 
*) S can also be written as S* = QQ J dx ^{tpip + 2W). 



Since the last two terms in Eq. (3-9) vanish under classical equations of motion, F = 
and {id — W")iIj = 0, relation (3-9) is consistent with the supersymmetry algebra, at least 
classically. 

After the Wick rotation, x -^ —ix and L -^ —L, we have the Euclidean action 

S = ldx!^^{d<pf + ^{Wr + Hd + W"),p-^F'^ (3-11) 

= -Q I dx]^{F + d<t) + W'), (3-12) 

and for the Hamiltonian H , 

\q {^{d(j) -W')]=H+\F{F + d(t)- W) - \^{d + W")^. (3-13) 

So far everything has been for the continuum. We now construct a lattice formulation of 
the above system on a finite-size lattice, 

yl = {xG aZ I <x </5}, (3-14) 

where a denotes the lattice spacing. First, we fix the lattice transcription of the time deriva- 
tive d. As a possible choice, we adopt the forward difference 

dj{x) = j{x^a)-j{x), (3-15) 

which does not lead to species doubling. The lattice counterparts of the Q transformation 
are then defined as 



[x)=i){x), QV^(x) = 0, (3-16) 

Qi){x) = F{x)-d4){x)-W\<l){x)), QF{x) = di){x)^W"{(\){x))i){x). (3-17) 

Finally, the lattice action is defined with an expression analogous to Eq. (3-12): 

S^-QY\ {^{^) (Fix) + d<Pix) + W^'(0(x)))} (3-18) 






^{W'{(j){x)))^ + ^{x) {d + W"{(j){x))) ^{x) 



-F{xY + W'{^{x))d^{x) \, (3-19) 



where 



+^{x = 0), for the periodic boundary condition, 
^(x = p) = { (3-20) 

-tjj{x = 0), for the antiperiodic boundary condition. 



9 



Note that all lattice variables are dimensionless. In the above lattice formulation with 
the periodic boundary condition, the Q-symmetry is manifestly preserved because the Q- 
transformations in Eqs. (3-16) and (3-17) are nilpotent, Q^ = 0. The invariance under Q 
is, however, broken because it is impossible to define a corresponding Q transformation 
on lattice variables such that the algebra {Q, Q} = —2d still holds. In fact, this lattice 
action is basically identical to the one described in Refs. 23), 27), 29), 35), 36), 45). (See also 
Ref. 30).) In some of these references, it has been shown that the Q-symmetry is restored 
in the continuum limit. 

As a Hamiltonian in this lattice formulation, following the discussion in the previous 
section and in view of Eq. (3-13), we use H(x) = iQQ{x)/2, where*^ 

Q{x) = --^{x) {id(l){x) - iW'{(/){x))) (3-21) 

is a lattice analogue of the Noether charge. The explicit form is 

H{x) = -^d(f){x)d(/){x) + ^{W'{(/){x))f - ^^{x) {d - W"{(l){x))) iIj{x) 

^(Jj ^(Jj ^Ui 

+ ^F{x) {d(l){x) - W'i(l){x))) . (3-22) 

The naive continuum limit of H{x) differs from the (imaginary-time) Hamiltonian in the 
continuum theory by terms vanishing under the classical equations of motion. **•* As discussed 
in the previous section, this lattice Hamiltonian has a correct zero-point energy in the sense 
that 

/ n rf0(x) dF{x) #(x) ckp{x) He~^ = 0, (3-23) 

which follows from the Q-exactness of H and the Q-invariance of the action S and of the 
integration measure. ***•* (Recall Eq. (2-9).) Thus this choice of the Hamiltonian is consistent 
with the topological nature of the Witten index, Eq. (2-3), with finite ultraviolet and infrared 
cutoffs. 

The numerical study of the present lattice model is not so difficult, because it is possible 
to obtain a closed expression of the fermion determinant in terms of (p. That is, 

det {-d - W"m = l[{l- W"{(l^{x))} T 1, (3-24) 

xGA 



*> We supplemented a factor of 1/a to adjust the physical mass dimension; recall that all lattice variables 
as well as the Q transformation are dimensionless. 

**' It is interesting to note that the expectation value (with the antiperiodic boundary condition) of this 
difference vanishes: (^^(F + dcp - W) - i^(9 + W")^;)^^^^ = 0. 

***■* Assuming that the integral /pg^ Hajeyi d<l){x) dF{x) d'ip{x) dip{x) Q e~^ is finite. This is certainly true 
if |M^'(±oo)| = +(X). 

10 



where the upper sign corresponds to the periodic boundary condition and the lower cor- 
responds to the antiperiodic boundary condition. Thus, after integrating over fermionic 
variables and the auxiliary variable F{x), we have the effective action 

Ses[<iy] = Yl {ld(j){x)d<l>{x) + ^-{W'{<j>{x))y + W'{<j>{x))d<j) 
l[{l-W"{<l^{x))}Tl 



xeA 



-In 



xgA 



(3-25) 



Note that the fermion determinant (3-24) is real for a real superpotential W but it is not 
necessarily positive definite. We thus must include the sign of the determinant, 

s[0] = sign (det {-d - W^"(0)}) = sign (nil- W"{(f){x))} ^ 1 j , (3-26) 

\xeA J 

as a reweighting factor in the functional integral. For example, the expectation value of the 
Hamiltonian H is given by the ratio of 

^f [YI d<P{x) Hs[(j)\ e-^-ffl'^l (3-27) 

to the partition function 

Z = AT /" JJ d(\){x) s[0] e~^'=«l'^'. (3-28) 

■^ xeA 

Recall that, when supersymmetry is spontaneously broken, the normalized expectation value 

with the periodic boundary condition cannot be defined because Z-p-QQ = 0. With the 

antiperiodic boundary condition, expectation values are always meaningful and the following 

substitutions can be made: 

(^(^)'>aPBC = -1' (^(^))aPBC = (3-29) 

for the auxiliary variable, and 

(^(x) {d + W"{<P{x))) ^(a;)>^p3c = -1 (3-30) 

and 

for fermionic variables. It is then straightforward to implement the hybrid Monte Carlo 
algorithm*^^) with the effective action (3-25) and compute (iJ(a;))aPBC- 

11 



Now, as a definite example in which supersymmetry is spontaneously broken, we con- 
sider*) 

1 1 

•^continuum = ^ ^^0continuum ' ofi'rcontinuum- [o-oZ) 

Since the parameter m has the mass dimension 1, we will measure all dimensionful quan- 
tities in units of m. For example, the lattice spacing is measured by the dimensionless 
combination am. If we introduce the dimensionless coupling constant as 

A = ^, (3-33) 

the superpotential in terms of the lattice variables reads 

W{(l){x)) = -{am)(P{xf + hamf/^\(p{xf. (3-34) 

Before proceeding to the numerical study, it is instructive to see how our method works 
for the free theory A = 0, for which supersymmetry is not spontaneously broken in the 
continuum theory and the lattice model (3-19) is solvable. From Eqs. (3-34) and (3-26), we 
see that the partition function Z (3-28) is nonzero for both boundary conditions (for a ^ 0). 
It is therefore meaningful to consider the expectation value for both boundary conditions 
and it is not difficult to see that 

{H{x))pBc = 0, (3-35) 

m(l — amy^"-~^ mil — amy ^""'^ ,^ ^„s 

^^(^^^^^^^ = 1 - (1 - amr/^ + 1 + il-amr/^- ^''^'^ 

In the continuum limit a -^ 0, the latter, in fact, reproduces the expectation value of the 
energy in this supersymmetric harmonic oscillator: 

hm(i/(x))aPBC = ^3^z^ + ^^^z^. (3-37) 

We thus have, in the large-time limit f3 ^ oo, 

Eo = lim lim(if(a:))aPBC = 0, (3-38) 

and infer that supersymmetry is not spontaneously broken. 

Now we turn to the Monte Carlo study of the model (3-34) with A 7^ (supersymmetry is 
dynamically broken in the target continuum theory). In the following results, we set A = 10 
and, for each set of parameters, we used 10^ statistically independent configurations. 



*^ The potential energy V{<l)) — W'icj))"^ /2 is a double-well type with two minima V" = at = and 
(j) — —m/g, and the height of the potential barrier is ?TT,*/(32g^) = m/(32A^). From this, the spontaneous 
supersymmetry breaking in the present system will be rather difficult to be observed numerically for weak 
couplings for which the supersymmetry breaking is caused mainly by quantum tunneling. 

12 
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Fig. 1. Histogram of the sign of the fermion determinant s[(j)] (3-26) (bold hne; left axis) and the 
Hamiltonian reweighted by the sign of the determinant Hs[(j)] (broken line; right axis) for the 
model (3-34) with the periodic boundary condition. A = 10. The lattice spacing is am = 0.1 
and the physical temporal size of the system is /3m = 1.6. The number of configurations is 10 . 



First, let us see the case of the periodic boundary condition. For this boundary condition, 
the sign of the fermion determinant (3-26) may change depending on the configuration, and 
in fact, as Fig. 1 shows, positive and negative fermion determinants appear at almost equal 
rates. This implies that the partition function (3-28), that is, the Witten index (2-2), is 
almost zero. This is perfectly in accord with the fact that supersymmetry is spontaneously 
broken in the target theory. In Fig. 1, we plotted also the distribution of the Hamilto- 
nian (3-22) reweighted by the sign of the determinant, Hs[(f)]. It spreads on negative as well 
as positive sides and the average (~ —0.003) is consistent with zero within statistical error 
(~ 0.02). This is again consistent with the fact that the average of the Hamiltonian (3-27) 
is merely the /9-derivative of the Witten index, Eq. (2-3). This assures us of the validity of 
our method because the construction of H ensures Eq. (2-3), as shown in Eq. (3-23). 

If we switch the boundary condition to antiperiodic, things drastically change. As Fig. 2 
shows, now the distribution of the sign is significantly asymmetric and the partition func- 
tion (3-28) becomes nonvanishing. This implies that we can give a definite meaning for the 
expectation value normalized by the partition function Z^p^q. In this way, we see numeri- 
cally that the effect of the boundary condition indeed survives even for large temporal size 
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Fig. 2. Histogram of the sign of the fermion determinant s[(f)] (3-26) (bold line; left axis) and the 
Hamiltonian reweighted by the sign of the determinant Hs[(p] (broken line; right axis) for the 
model (3-34) with the antiperiodic boundary condition. A = 10. The lattice spacing is am = 0.1 
and the physical temporal size of the system is I3m = 1.6. The number of configurations is 10*^. 

{j3m = 1.6 is actually a large size with the present value of the coupling constant; see below) 
when supersymmetry is spontaneously broken. 

With the antiperiodic boundary condition, we then measure the expectation value of the 
Hamiltonian as a function of the temporal size of the system (3. For various values of /3m, 
we measured {H {x)) i^p-qq / m for lattice spacings am = 0.1, 0.05 and 0.02. The number of 
configurations is 10^ for each set of parameters. Then, as shown in Fig. 3, we extrapolate 
{H (yx)) i^p-QQ / m to the continuum a = by a linear x^-fit.*^**^ In Fig. 4, we plot the continuum 
limit of the expectation value \mva-tQ{H{x))sP^c/m as a function of the physical temporal 
size of the system (3m. For /9m > 1^***) we have \m\a-^Q{H [x)) ^pbc — 1-1^^ and, from this. 



*' The statistical errors in Fig. 3 are one standard deviation. The errors in the linear x^-extrapolation 
were estimated from the range of fitting parameters that corresponds to a unit variation of x^ . These remarks 
also apply to the results in Fig. 5. 

**' Wc found that a quadratic function of the form a{amy' + /3, which gives a somewhat larger 
{H (x)) aPBC / ^TT' at a = 0, provides a better fit. Although this form of the fit function might be suggested 
theoretically (i.e., the residual lattice artifact is 0{a^) instead of 0{a)), we stick to a simple linear fit to 
avoid a possible criticism that nonzero Eq is an artificial consequence of the fit. 

***) por strong couplings A 3> 1, it is easy to see that energy eigenvalues of the present system scale 
as X^/^m. Thus the difference between the first excited state and the ground-state would be X^/^m times 
a number of 0(1). This observation suggests that the expectation value of the Hamiltonian (with the 
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0.12 



Fig. 3. Linear extrapolations of {H {x)) sP'bc / "^ to the continuum a = for various values of [5m. 
A = 10. The errors are only statistical ones. 




Fig. 4. The continuum limit of the expectation value of the Hamiltonian, \vca.a^o{H{x))sp-Qc/'>^-, 
as a function of the physical temporal size of the system (5m. A = 10. The errors are only 
statistical ones. We have also plotted the exact ground-state energy Eq/tji = 1.27616 and the 
analytic expression for the A = case, Eq. (3-37). 
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we infer that supersymmetry is spontaneously broken. *•* Indeed, this is the correct answer. 



3.5 



2.5 



1.5 - 

1 
0.5 


-0.5 



"1 1 1 

aPBC I — H- 

PBC ^--x- 



^-?-^-|-*l- 



«-■ 



1 



|3m 



Fig. 5. Continuum limit of the expectation values of the Hamiltonian, lima^o(-f^(3^))aPBc/'^ ^-nd 
lima^o{H{x))pBc/''TT-: ^s a function of the physical temporal size of the system p-m. The errors 
are only statistical ones. 

Next, as an interacting case in which supersymmetry is not spontaneously broken, we 
consider 



W, 



continuum 



Tfi (h 

r continuum' 



(3-39) 



In this case, we observed that s[0] (3-26) has almost always a definite sign for both boundary 
conditions and thus both {H{x))pbc and {H{x))aPBC can be considered. In Fig. 5, we plot 
the continuum limit of these quantities as a function of Pm. The results are obtained by 
extrapolation to the continuum a = with a linear x^"fif (like Fig. 3) of data computed 
at am = 0.1 and 0.05. The number of configurations is 10^ for each set of parameters. 
The figure shows that, in this case for which supersymmetry is not spontaneously broken, 
liuia-^o{H{x))pBc is consistent with zero for all temporal sizes (recall Eq. (2-4); in the present 
model, Tr(— 1)-^ = 1) and lima^o(-^(2;))aPBC approaches zero as the temporal size of the 



antiperiodic boundary condition) exponentially approaches the asymptotic value at (3m — oo around f3m > 
A-2/3 ~ 0.2 for A = 10. 

*^ In Fig. 4, we plotted also the exact ground-state energy Eo/m — 1.27616 that was obtained by 
numerically diagonalizing the corresponding Hamiltonian operator (we used the method of Ref. 47)). The 
discrepancy of our Monte Carlo results for (3m > 1 with this exact result can be understood as a systematic 
error associated with a linear extrapolation to the continuum limit. 
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system is increased. From this, we conclude that -Eq = is within the error. 

In summary, we have observed that our method works perfectly well in the present su- 
persymmetric quantum mechanics. One can certainly observe whether or not the dynamical 
supersymmetry breaking takes place by our method. 

§4. Two-dimensional J\f = (2,2) super Yang-Mills theory 

The two-dimensional A/" = (2, 2) super Yang-Mills theory is obtained by a dimensional 
reduction of the four-dimensional J\f = 1 super Yang-Mills theory. *) This seemingly simple 
supersymmetric system, however, defies a straightforward low-energy description for several 
reasons. First, two global U{1) symmetries in this system cannot be spontaneously broken 
in two dimensions and a description by using the Nambu-Goldstone fields is impossible.**^ 
Second, there is no controllable parameter, other than the number of colors N^ of the gauge 
group SU{Nc). (The two-dimensional gauge coupling g simply provides a mass scale, just 
like ylqcD-) The 1/Nc expansion is nontrivial because the gaugino and scalars belong to 
the adjoint representation. Finally, the classical potential energy of scalar fields possesses 
noncompact flat directions and there are an infinite number of degenerated classical vacua. 
This classical degeneracy is not lifted upon quantum corrections to all orders of perturbation 
theory. 

In our present context, the last point above (noncompact flat directions in the classical 
potential) is an obstruction to the determination of the Witten index. In the weak coupling 
approximation, zero-momentum modes without potential (constant degrees of freedom along 
flat directions) produce a continuous spectrum starting at zero. This makes the counting 
of zero-energy states, and thus the determination of the Witten index in the weak cou- 
pling approximation, awkward. A similar situation arises in the three-dimensional A/" = 2 
super Yang-Mills theory (that can also be obtained by a dimensional reduction of the four- 
dimensional A/" = 1 super Yang-Mills theory). However, in this three-dimensional model, if 
the gauge group is SU{Nc), one eliminate zero-momentum bosonic modes by imposing the 
twisted boundary conditions for two spatial directions and obtain Tr(— 1)^ = 1.49)***) 'jj-^jg 
trick of the twisted boundary conditions, unfortunately, does not work in two dimensions. 
The correct value of the Witten index, or even whether it is well deflned or not, is therefore 



*) In what follows, we assume that the gauge group is SU{Nc). 

**) Nevertheless, it is possible to show that a correlation function of Noether currents associated with the 
C/(l) symmetries possesses a massless pole, to all orders of perturbation theory.'*®^ 

***' Incidentally, this is a good example of the Witten index generally not being preserved under dimen- 
sional reduction, because Tr(— 1)^ = Nc for the four-dimensional A/" = 1 super Yang-Mills theory. See also 
§ 4of Ref. 2). 
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unknown for the two-dimensional M = (2, 2) super Yang-Mills theory. It is consequently not 
known whether supersymmetry is dynamically broken in this system. 

Under this situation, Hori and Tong^-* conjectured that dynamical supersymmetry break- 
ing occurs in this system, on the basis of the counting of the number of ground-states in 
the two-dimensional A/" = (2,2) supersymmetric gauge theory with fundamental chiral mul- 
tiplets, combined with a decoupling argument. In what follows, we numerically investigate 
this possibility of dynamical supersymmetry breaking by directly measuring the ground-state 
energy density in Euclidean lattice gauge theory. 

4.1. Hamiltonian density in the continuum theory 

The Lagrangian density of the two-dimensional M = (2, 2) super Yang-Mills theory in 
the Minkowski spacetime, in terms of the twisted basis of spinors,^'^-''^^-' is given by 



+ 7V[<P, V] + X[0, X] - V^m[0, ^^^] + '^ixiiDo^pi + Ditpo) - tpoDoT] - iipiDiT] 



1 

(4-1) 



where all fields are SU{Nc) Lie algebra valued and scalar fields and are combinations of 
two real scalar fields, = X2 + iX^ and = X2 — iX^, respectively. Fqi = OqAi — diAo + 
i[Ao,Ai] is the field strength in two dimensions. The covariant derivatives D^ are defined 
with respect to the adjoint representation D^ip = df^ip + i[A^, ip] for any field (/?. The index fi 
runs over and 1. Note that, in the above convention, the bosonic fields A^, and have 
the mass dimension 1 and the fermionic fields -0^, x ^''^^ V have the mass dimension 3/2, 
because the gauge coupling constant in two dimensions g has the mass dimension 1. 

The action J (Px C with the periodic boundary condition is invariant under four su- 
persymmetry transformations. Among them, what is relevant to us is Q and Qo- The 
Q-transformation is given by 



QAq = iipo, 




Qipo = ^00, 


QAi = V-i, 




Q^i = iD^<p 


Q0 = O, 






Qx = H, 




QH =[(!>, xl 


Q4> = ri, 




Qv = [0, 0! 


and Qq is (see, for example, Ref. 52)) 








Qov = 


-2Do0, 



(4-2) 



QqAi = -X, QoX = iDi(j), 

Qo0 = 0, 

QoV^i = H- 2Foi, QoH = -[0, ^i] - 2DoX - ^DlV■ 

Qo0 = -2^o, Qo^o = ^[0,0]. (4-3) 

One then finds that these transformations satisfy 

Q' = 5^, Ql = -S^, {Q,Qo} = -2do-2t5Ao, (44) 

where (5^ denotes the infinitesimal gauge transformation with the parameter ip. These differ 
from the off-shell supersymmetry algebra in the twisted basis, Q^ = Ql = and {Q, Qq} = 
— 2(9o, by gauge transformations because we are working with the Wess-Zumino gauge. 

A crucial property of this system, which allows a simple lattice formulation, is that the 
action S = f (Px C is Q-exact.^^^'^^^ 

S = q\J d'^x tr |-^r7[0, 0] + 2xFoi - xH + ^oDo^ + ^Vi^i^} (4-5) 

In this form, with the relations (4-4), the invariance of the action under Q and Qq transfor- 
mations is easily seen.*) 

Now, from the Lagrangian density (4-1), we obtain the Hamiltonian density 7i by the 
Legendre transformation. After (trivially) eliminating redundant fields by using second-class 
constraints, 7i is given, in terms of fields in the Lagrangian, by 

1 f 1 . , -V 



n = -tT<>- [0, 0]2 + H^ + Dqcj^DqcI) + Di0Di0 

- ^^[0, V] - X[0, X] + i'i^\4>, V^m] - 2ix^i^o + ii^iDiT] 
-2tr{Ao6^}, (4-6) 

up to a spatial total derivative, where Q is the Gauss-law constraint: 

^ = 1 |DiiJ + ^[0,Do0] + ^[0,/^o0] +H^1,X} + ^K^o}} . (4-7) 

From the off-shell supersymmetry algebra {Q, Qo] = —2dQ, one might expect that the rela- 
tion QJ'q/2 = Ti holds, where Jq is the time component of the Noether current associated 
with the Qo-syinmetry 

Jo = \ tr |ir/[0, 0] + 2xH + 2i,QD^ - 2ii,rD;A . (4-8) 



*) Note that S can further be written a& S — QQo-^ I '^^^ tr {~^(j)Do(j) — V'lx}- 
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In reality, 



\qJ^ = n + 2tT {AoG} + ^ tr {^0 (-2[0, H + 2^^iX - ^o^) } (4-9) 

^ 9 



up to a spatial total derivative. Compare this with Eq. (2-8). In Eq. (4-9), the last two 
terms are proportional to classical equations of motion and they can be expressed as Aq-^S 
and 4'ojj;s:S.*^ The expectation value of these two expressions may be obtained (after gauge 
fixing) as a Jacobian associated with the transformations Aq -^ AQ + aAp, and ipQ -^ ipQ + aipQ, 
respectively. Such a Jacobian is generally ultraviolet divergent.**^ 

Thus, in view of Eq. (4-9) and following our general prescription, we adopt 7i = QJq/2 
as the Hamiltonian density in our functional integral formulation. As can be seen from 
Eqs. (4-9) and (4-6), this Ti is, moreover, gauge invariant. 

Thus, going to the Euchdean space by Xq -^ —ixo, Aq -^ iA^, Dq -^ iDq and C -^ — £, 
we have the Euclidean action 

S = Q^j(fx tr |ir/[0, 0] - ix<^ + XH- i^^^D^A , (4-10) 

where ^ = 2Foi and 

QA^ = ^^, Qi)^ = iDf^cp, 

QX = H, QH=[<P,x], 

Q4> = v, Q^=[0,0], (4-11) 

and the Hamiltonian density 

1 fi 

These are the basic relations for our Euclidean lattice formulation. 



'H = Q-ii{ -77[0, <P]+xH + i^oDo(P - i^iD^cj) } . (4-12) 



4.2. Manifestly Q-invariant lattice formulation 

This subsection is a brief summary of a lattice formulation of the two-dimensional M = 
(2,2) super Yang-Mills theory proposed in Ref. 13). For full details, we refer the reader to 
Ref. 13). We consider a two-dimensional rectangular lattice of the physical size j3 x L: 

A= [xeaZ^ \Q<xq< (3, < xi < L} , (4-13) 



*•* We define color components of fields as ip = X]a=i ^"'T"-, where T" are generators of SU{Nc). 
**' With a lattice regularization, for example, the one described in the next subsection, the expectation 
value of the latter is —i{N^ — l)/^^- The expectation value of the former is +i{N^ — l)/a^ after gauge fixing 
and thus, quite interestingly, the expectation value of the last two terms of Eq. (4-9) vanishes. 
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where a denotes the lattice spacing. All fields except the gauge potentials are put on sites 
and the gauge field is expressed by the compact link variables f/(x, /i) G SU{Nc). 
As a lattice transcription of the Q-transformation (4-11), we define 

QU{x,j2) = iiljf,{x)U{x,j2), 

Q^J^ix) = iij^{x)ij^{x) - i {(f){x) - U{x, fi)(f){x + afi)U{x, /i)"^) , 

g0(a;) = 0, 

Qx{x) = H{x), QH{x) = [<p{x), x{x)], 

Q4>{x) = r]{x), Qr]{x) = [(/){x) ,^{x)] (4-14) 

{fi implies a unit vector in the /i-direction). It can be confirmed that Q^ = S^, where 6^ is 
an infinitesimal gauge transformation on the lattice with the parameter (f){x). The lattice 
action is then defined by an expression analogous to Eq. (4-10): 

S = Qa^\2 (o,{x) + 02ix) + Os{x) + ^tT{x{x)Hix)}) , (4-15) 



xGA 



where 



0^{x) = -^tTl^r]{xMx),(l^{x)]\, (4-16) 

^2(a:) = ^tT[-ix{x)Six)} , (4-17) 

^3{x) = -^-^ tr J i V ^^(x) (0(a;) - U{x, /i)0(x + afi)U{x, fi)-^) } . (4-18) 

In Eq. (4-17), ${x) is a lattice counterpart of the field strength and is defined from the 
plaquette variables 

U{x, 0, 1) = U{x, 0)U{x + aO, l)U{x + al, 0)-^U{x, l)'^ (4-19) 



as 



S(x) = '^ ^, ^x) = -i\U{x,0,l)-U{x,0,l)-^], (4-20) 

^ ' l-i||l-t/(x, 0,1)11' ^ ^ L V > > ; V , , ; J, v ; 



where the matrix norm is 

1/2 



= [ti {AA^JY^' (4-21) 
and the constant e is chosen in the range 

< e < 2V2, for Nc = 2, 3, 4, (4-22) 

0<e<2v/iVeSin(^J , iorN^>5. (4-23) 
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From the Q-exact form (4-15) and the nilpotency of Q, Q^ = S,/,, the lattice action is 
manifestly invariant under the Q-transformation (4-14).*^ 
After the operation of Q, the lattice action becomes 
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If 1 -^ 

Bi[x} + > , £Fi[x) + — -- tr <^ H{x) - -^^tlW 



1=1 



(4-24) 



where we have noted that only the traceless part of ^(x) 



<pTL{x)=^{x)-—tr{<P{x)\l, 



(4-25) 



appears in the action, because the auxiliary field H{x) is traceless. Each term of the action 
density is given by 



Cn^(x) = — r^tr < -[(f){x) , (p^x)]"^ 



4^2 --141 






a 9 
1 



1 
A" 



tr<! -<pTL(a; 



4^2 - I 4 



tr<^ (0(x) - U{x, O)0(x + aO)?7(x, 0)-^) 



(4-26) 
(4-27) 



X {(j){x) - U{x, O)0(x + aO)f/(x, 0)"^) 



(4-28) 



CnA(x) 



a^g^ 



tr<^ (0(a;) - U{x, l)0(x + ai)f/(x, 1)"^) 



X {(pix) - U{x, l)(j){x + ai)U{x, ly^) 



(4-29) 



and 



Cfi{x 

Cf2{x 
Cf3{x 
Cf4^{x 

Cpsix 



-— tr I --r]{x) [(t){x),r]{x)] \ , 

1 

tr{-x(a;)[0(a;),x(a;)]} 



aW 



a^g'^ 



a^g'^ 



a^g^ 



ii {-i)Q{x)i)Q{x) [(t){x) + U{x, O)0(x + aO)f/(a;, 0)"^) } , 
tr {-?/'i(a;)V^i(a;) (0(a;) + f/(a;, l)0(a; + al)U{x, 1)"^) } , 
tr|ix(a;)<3^(a;)}, 



(4-30) 
(4-31) 
(4-32) 
(4-33) 
(4-34) 



*' Another interesting property of the present lattice formulation is that one global U{\)b. symmetry is 
manifestly preserved. ^^^ 
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1 
Cpeix) = -r^ tr {-#o(a;) (vix) - U(x, 0)r](x + aO)U(x, 0)"^) } . (4-35) 

'^Frix) = -^ tT{-iiJi{x) (r]{x) - U{x, l)r]{x + ai)U{x, 1)-') } . (4-36) 

Note that all lattice fields in the above expressions are dimensionless. For all fields other 
than fermionic fields, periodic boundary conditions on A are assumed. For fermionic fields, 
ip = [iIjh,XjV)j depending on whether the temporal boundary condition is periodic (PBC) 
or ant iper iodic (aPBC), we set 

{+ip{xo = 0,a;i) for the periodic boundary condition, 
(4-37) 
—■ip{xo = 0,Xi) for the antiperiodic boundary condition, 

while the spatial boundary condition is always taken to be periodic. 
With the lattice action (4-24), the partition function is defined by 



Z = X f dfxe-'^, (4-38) 



where the integration measure is defined (writing (f){x) = X2{x)+iX3{x) and 0(a;) = X2{x) — 
iXsix)) by 

^/^ - n ( n ^^(^' *") ) fl dX^(^) dX^{x) dH^{x) I Yl dij^ix) j dx^'ix) dri^ix) (4-39) 

X&A \^l=0 J a=\ \M=0 / 

in terms of color components of fields. dU{x, fi) is the standard Haar measure. Note that 
the integration over the auxiliary field H{x) can readily be performed because it is gaussian. 
The invariance of this measure under the Q-transformation is noted in Ref. 17). 

The denominator in Eq. (4-20) needs an explanation. Without that factor, the lattice 
action for the gauge field is the "double-winding plaquette type'"^^^ and the action possesses 
many degenerate minima which have no continuum counterpart. Because of the denominator 
of Eq. (4-20), the action (4-24) diverges as ||1 — U{x, 0, 1)|| — ^ e at a certain site x. Precisely 
speaking, the above construction of the action is applied only for configurations with 

||l-f/(a;,0, 1)11 <e, for Vx G yl, (4-40) 

and, otherwise, i.e., if there exists x E A such that ||1 — U{x, 0, 1)|| > e, we set 

S = +oo. (4-41) 

In this way, the domain of functional integral (4-38) is effectively restricted to configurations 
specified by Eq. (4-40). It can then be shown that U{x,fi) = 1 is (up to gauge transforma- 
tions) a unique minimum of the action within the integration domain. This procedure to 
solve the problem of degenerate minima moreover does not break the Q-symmetry.^^-' 
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With the above construction (and with the periodic boundary condition), one fermionic 
symmetry Q is manifestly preserved on the lattice. The price to pay is that the Pfaffian of 
the Dirac operator, resulting from the integration over fermionic fields, is generally complex. 
Since the corresponding PfafSan in the target continuum theory is real and positive semi- 
definite, the complex phase must be a lattice artifact. That is, we expect that the imaginary 
part of the lattice Pfaffian diminishes as we approach the continuum limit. We will later 
confirm this expectation numerically. 

In Refs. 12) and 13), the restoration of the invariance under a full set of supersymmetry 
transformation in the continuum limit has been argued on the basis of perturbative power 
counting. It is certainly desirable to confirm, however, this restoration nonperturbatively 
by observing the supersymmetric Ward-Takahashi identities. More definitely, one should 
examine the total divergence of two-point (denominator-free) correlation functions containing 
the supercurrent; so far this analysis has not yet been carried out. 

4.3. Lattice transcription of Hamiltonian density 

Now, as the definition of a Hamiltonian density on the lattice, we follow the prescription 
n{x) = QJ^{x)/2 suggested by Eq. (4-12), where 

Joi^) = -^tT\lr]{xMx),^ix)]+2xix)Hix) 
a, g [^ A 

- 2iiJo{x) (0(a;) - U{x, O)0(x + a6)U{x, 0)"^) 

+ 2#i(a;) {^{x) - U{x, l)0(a; + ai)U{x, 1)"^) I (4-42) 



is a lattice transcription of the Noether current J'q in Eq. (4-12). The explicit form of the 
Hamiltonian density is then given by 

H{x) = Cbi{x) - Cb3{x) + Cb4{x) 

+ Cfi{x) + CF2ix) - CF3ix) + Cf4{x) - Lf%{x) ^ Cfy{x) 

+ ^tT{Hix)}\ (4-43) 

From the Q-invariance of the lattice action and of the integration measure, we thus have 

f dfiHix) e''"^ = f dfxQ(-J^{x)e-A =0, (4-44) 

ipBC JPBC \2 / 

assuming that the integral Jp^^dfiJ'Q{x)e~^ is finite, and this reproduces the topological 
property of the Witten index, Eq. (2-3). As already discussed, we regard this property as a 
guiding principle for choosing the origin of the energy (density). 
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We will measure the ground-state (vacuum) energy density Sq by 

lim \im{n{x)),PBC = lim lim -^aPBC ^/^^(^) ^ ^ g^ (4.45) 

and judge that dynamical supersymmetry breaking occurs if £^o > and but not if £q = 0. 
For Eq. (4-43), an integration over the auxiliary field H{x) can be performed to obtain 

\ " ' / aPBC 

= ^(A^c-l)^-(^B2(x)),pBc- (4-46) 

Thus, in actual numerical simulations, we can (suppressing the subscript aPBC or PBC) use 

{n{x)) = {Cbi{x)) - {Cb2{x)) - {Cb3{x)) + {Cb4{x)) 

+ {Cfi{x)) + {Cf2{x)) - {Cf3{x)) + {Cpiix)) - {Cfq{x)) + {Cprix)) 

+ l(Nc-l)^r (4-47) 

We can argue that the expectation value (4-47) with aPBC is ultraviolet finite and pos- 
sesses a well defined continuum limit to all orders of (lattice) perturbation theory. In the 
present super-renormalizable model, a possible ultraviolet divergence in Eq. (4-47) arises 
from either (i) one-loop two-point functions of bosons contained as sub-diagrams or (ii) one- 
loop diagrams that are formed by a self-contraction of kinetic terms. The former (i) case 
is potentially logarithmically divergent but the divergence is cancelled among boson loops 
and fermion loops. For PBC, this cancellation can be shown by using the Q-symmetry. 
For aPBC, the cancellation still holds because the change in the boundary condition does 
not influence the coefficients of the logarithmic divergent pieces. In the latter case (ii), the 
divergence is quadratic and it is common to the free theory. We may thus examine the 
expectation value in the free theory and find that, even with aPBC, it is ultraviolet finite 
for any f3 > 0. Our numerical study described below in fact indicates that the continuum 
limit of the expectation value (7i(a;))aPBc is well defined. 

4.4. Monte Carlo study 

We numerically studied only the SU{2) gauge group. Our algorithm and the computation 
code, which was developed using FermiQCD/MDP,^^^'^^^ are almost identical to those in 
Ref. 20). We use the hybrid Monte Carlo algorithm*^^) to generate configurations in the 
quenched approximation. The effect of dynamical fermions is later taken into account by 
reweighting configurations by the Pfaffian of the Dirac operator. We do not introduce any 
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mass terms of fermions or bosons that would explicitly break the Q-symmetry. Although 
this is certainly a brute force method compared with a standard pseudo-fermion algorithm, 
its implementation is much simpler and the validity has been confirmed for one-point Ward- 
Takahashi identities.^^-* 

A direct calculation of the Pfafiian is very time-consuming. *) Thus we instead use a square 
root of the determinant that is obtained by LU decomposition. Expressing the determinant 
of the Dirac operator D in the form 

det{D} = re'\ -tt < ^ < vr, (4-48) 

we evaluate the Pfafiian by 

Vi{D} = ^e'^l^, (4-49) 

because (Pf{-D})^ = det{i5}. This prescription, however, reproduces a correct Pfafiian if 
and only if — 7r/2 < Arg(Pf{Z)}) < 7r/2. It is expected that this inequality is fulfilled in the 
continuum limit, because the Pfafiian in the continuum target theory is real and positive 
semi-definite. A direct calculation of the Pfafiian over a subset of our configurations (Fig. 6) 
clearly supports this expectation and justifies the above prescription. Note that our present 
lattices are much finer, compared with that in Ref. 20) where ag > 0.5. 

We stored statistically independent configurations for parameters summarized in Table I, 
where Nt and Ns are the number of lattice points for the temporal and spatial directions, 
respectively. The physical size of the spatial direction is fixed to be Lg = v2-** The 
parameter e in Eq. (4-22) is taken to be e = 2.6. We used the cold start and set all scalar 
fields to be zero at the initial configuration. As the initial thermalization, we discarded 
the first 10^ trajectories and then stored configurations at every 10^ trajectories (the auto- 
correlation time was 10-20 trajectories). 

To give an idea of the quality of our numerical simulation and to illustrate that the 
quantum effect of fermions is really taken into account, in Fig. 7, we plot the real part of the 
expectation value of the action density in Eq. (4-15) with the periodic boundary condition as 
a function of the lattice spacing ag. Since the lattice action density is Q-exact, its expectation 
value under the periodic boundary condition should be zero for any lattice spacing. The plot 
is certainly consistent with this. On the other hand, the expectation values in the quenched 
approximation are definitely not consistent with zero as they do not contain the effect of 
dynamical fermions. See also Fig. 2 of Ref. 20). 



*' It can be seen that the algorithm for the Pfaffian (appearing, for example, in Ref. 55)) is an 0{n'^) 
process for a 2n x 2n matrix, while the LU decomposition has an 0(n'')-process algorithm. 

**' Note that supersymmetry is not broken in the infinite volume if it is not with finite volume. "^^ This 
fact would justify our study with finite physical volume. 
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Fig. 6. Histogram of the complex phase of the Pfaffian Arg(Pf{D}) in radians obtained by direct 
calculation of the Pfaffian for sampled configurations. The plots are for two entries in the 
Nt/Ns = 1 column in Table I and the number of sampled configurations is 160 for both cases. 
The boundary condition is antiperiodic. 



Table 1. Number of statistically independent configurations we used for the cases with the an- 
tiperiodic boundary condition. Nt x Ns = 3 x 6 is the minimal-size lattice and 36 x 12 is the 
maximal. The physical spatial size is held fixed at Lg = \/2 = 1.4142. 
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Now, our main result in this paper is illustrated in Fig. 8. We plotted the continuum 
limit of the real part of the expectation value of the Hamiltonian density (4-43) with the an- 
tiperiodic boundary condition, lima^oR-e('^(2;))aPBC; as a function of the physical temporal 
size of the system f3g. For each /Sg, the continuum limit was obtained by a linear x^"fit, as 
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Fig. 7. Real part of expectation values of the action density (over g^) with the periodic boundary 
condition. The parameters are identical to those of entries in the Nx/Ns = 1 column in Table I, 
except that the number of configurations is 9,900 for each case. 
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lim^^o Re('H(x))aPBc/5^ and lim.a^oRe{7{{x))pBc/9'^, as functions of the temporal size f3g. 
The errors are only statistical ones. For the periodic boundary condition, the number of 
configurations is 9,900 for all cases. 
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depicted in Fig. 9.*-* For Pg ^ I, the expectation value rapidly approaches the asymptotic 
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Fig. 9. Linear extrapolations of Re('H(x))aPBc/s'^ to the continuum a = for various values of 
Pg. The errors are only statistical ones. 

value, that is, So/g"^, according to Eq. (4-45). We may estimate the asymptotic value in 
f3g ^ oo hj x^-fit using a constant. The use of four data points in Fig. 8 at (3g > 2 gives 
^o/g"^ = —2.0 ±1.7 and three points at Pg > 2.5 gives So/g^ = —3.0 ± 2.2. If we use an 
exponential function Aexp{—BPg) + C and all points in Fig. 8, we have So/g"^ = —2.2 ± 1.4. 
All these results on So/g"^ are consistent and, at least within one standard deviation, we do 
not observe positive vacuum energy density. We regard this as an indication of the fact that 
supersymmetry is not dynamically broken in this system. Of course, errors in our present 
result are large and we cannot exclude the possibility of supersymmetry breaking of 0(1) in 
So/g'^. Further reduction of statistical errors will allow us to conclude whether the scale of 
dynamical supersymmetry breaking is 0(1) or not.**) 

In the present lattice model with the periodic boundary condition, we observed that the 
Pfafiian of the Dirac operator is almost real positive (recall Fig. 6) and this implies that 

*■* The statistical errors in Fig. 9 are one standard deviation, obtained by jackknife analysis. Jackknife 
analysis is necessary because we are using the reweighting method, as explained in Ref. 20). The errors in 
the linear ^^-extrapolation are estimated from the range of fitting parameters that corresponds to a unit 
variation of x^ . 

**' For the present lattice model, we are currently developing a simulation code with the pseudo-fermion 
and the RHMC algorithm. ^^^ We hope this will enable us to reduce the statistical errors without substantially 
increasing the number of configurations. 



29 



-2pBc 7^ with finite lattice spacings (unlike the case of Fig. l).*-* We can thus consider the 
expectation values with the periodic boundary condition. In Fig. 8, we have also plotted 
the real part of the expectation values of the Hamiltonian density for various temporal sizes 
with the periodic boundary condition. Since the Hamiltonian density is Q-exact, all the 
expectation values with the periodic boundary condition must be zero, if one can define them. 
The plot is clearly consistent with this. This also supports the idea that supersymmetry is 
not broken in this system. If supersymmetry is spontaneously broken, the expectation value 
(7Y(x))pBc must be indefinite, as Eq. (2-4) shows. In our simulation, the expectation value 
is computed as {H^x) Pf{-D})qucnchcd/(Pf{-D})qucnchcd!^°'' and it can be indefinite only when 
(Pf{D}) quenched = in the coutinuum limit. We did not see such a tendency and obtained 
the plot shown in Fig. 8. Also, it is interesting to note that the overall feature of Fig. 8 is 
quite similar to that of Fig. 5, rather than that of Fig. 4. 

§5. Discussion 

The most direct way to observe the spontaneous supersymmetry breaking would be to 
examine the degeneracy of boson and fermion mass spectra through two-point correlation 
functions. Although this method is conceptually clear, a reliable exponential fit of two- 
point functions would require a rather large lattice extent. The method we propose in this 
paper is computationally much easier because it is based on the measurement of one-point 
functions, the expectation values of a Hamiltonian (density). A weakness is the ambiguity 
in the choice of the Hamiltonian in the Euclidean lattice formulation. In this paper, we gave 
a justification (on the basis of a topological property of the Witten index) of the choice, for 
lattice formulations that possess a manifestly preserved fermionic symmetry Q. In any case, 
this is the first work of a direct investigation of the spontaneous supersymmetry breaking in a 
gauge field model (for which the Witten index is unknown) by numerical simulation. Before 
the recent developments in the lattice formulation of supersymmetric gauge theories, ^-'"^•* one 
could not even imagine such a study feasible. 

One may ask the extent of the applicability of our method. We already have a lattice 
formulation^^^ of the two-dimensional J\f = (4, 4) super Yang-Mills theory in which two 



*' In the quantum mechanical system (3-32) in which supersymmetry is spontaneously broken, the Witten 
index ZpBC (2-2) becomes zero because the fermion determinant is not positive-definite, as shown in Fig. 1. 
In our present lattice model for the two-dimensional M — (2, 2) super Yang-Mills theory, the PfafBan is 
almost real positive, and thus, -ZpBC 7^ with finite lattice spacings. One might then think that this latter 
fact alone is sufficient to conclude that supersymmetry is not spontaneously broken in this system. Although 
this argument is not quite correct, because there is a possibility that the coefScient A/pBC, and thus ZpBc, in 
Eq. (2-2) becomes zero in the continuum limit, it certainly indicates that dynamical supersymmetry breaking 
is unlikely in this system. 
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fermionic symmetries are exactly preserved. Thus, it should be possible to study possible 
dynamical supersymmetry breaking in this theory in a similar manner. For other supersym- 
metric field theories (except Wess-Zumino-type models), strictly speaking, we do not have a 
lattice formulation with manifestly preserved fermionic symmetry. Further study is needed 
on the lattice formulation of these theories (including physically interesting models, such as 
two-dimensional supersymmetric nonlinear sigma models, three-dimensional supersymmet- 
ric pure Yang-Mills theories, and supersymmetric gauge theories with matter multiplets). 
For related works, see Refs. 57)-60). 

Suppose that we have a lattice formulation in which the lattice action S and the in- 
tegration measure dfi are manifestly invariant under a fermionic transformation Q. It is 
quite conceivable that, for typical models in dimensions higher than two, this manifest Q- 
invariance alone is not sufficient to ensure automatic restoration of the invariance under a 
full set of supersymmetry transformations. One would then have to supplement a counter 
term AS to the original lattice action. However, it is also conceivable that AS is invariant 
under Q, because the original lattice regularization preserves a manifest Q-invariance. If 
this is true, it is again natural to adopt the prescription for the Hamiltonian H = iQQ/2 
because the relation 

df^He-'-'^'= [ df^'-QQe-'-^' = [ df,Q ('-Qe''-^'] = 0, (5-1) 

PBC JPBC ^ J PEG \^ / 

which corresponds to the topological property of the Witten index, still holds. 

Can we not do anything if the Q-invariance is not manifest in the lattice formulation 
that is adopted? A natural idea is to take an arbitrarily chosen Hamiltonian H and then 
subtract a constant from it, H = H — c, such that relation (2-3) holds. It is easy to see that 
this requirement implies*-* 

H = H- -^P^^ ^ _^ . (5-2) 

JpBc ^^ ^ 

However, when supersymmetry is spontaneously broken, the denominator of the second term 
of Eq. (5-2) would vanish (because it is proportional to the Witten index) and thus, unfortu- 
nately, it appears that formula (5-2) itself cannot be used for cases in which supersymmetry 
is spontaneously broken. We certainly need a more elaborate idea. 

*' Incidentally, this formula reproduces a prescription for the origin of the energy in Refs. 31) and 
32) in which the thermal average of the energy in one-dimensional supersymmetric Yang-Mills theories is 
numerically studied. The thermal average of the energy is given by (the minus) the /3-derivative of the 
thermal partition function (2-5). In one-dimensional supersymmetric Yang-Mills theories, by rescaling the 
imaginary time and dynamical variables, one sees that a Hamiltonian H is simply given by —3/(3 times the 
Euclidean action up to an additive constant. If one substitutes this H into Eq. (5-2), one ends up with 
the formulas in Refs. 31) and 32). Assuming that Jp^„diie^^ ^ 0, the second term in Eq. (5-2) can be 
evaluated by the lowest-order perturbation theory. 
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Our main aim in this work was to examine a possible spontaneous supersymmetry break- 
ing from the ground-state (vacuum) energy obtained by hm/3^oo(-f^)aPBC- It is nevertheless 
important to study {H)aPBC with finite (3 because it contains useful information on the 
energy spectrum of excited states. That is, when there is an energy gap AE between the 
first excited state and the ground-state, the decay of {H),^^^q for /? -^ oo is exponential, 
~ exp(— /3Z\i5), whereas when the spectrum is continuous starting at zero and the density 
of states behaves as p{E) ~ E^^^, the decay of (i/)aPBC for /3 — * oo is power-like, ~ z///9. 

The behavior in Figs. 4 and 5 appears to be consistent with the exponential decay ex- 
pected for quantum mechanical systems with discrete spectra. For the two-dimensional 
M = (2, 2) super Yang-Mills theory in Fig. 8, is the decay an exponential or power one? The 
error bars in the figure are too large for a reliable fit and we reserve this study for a future 
work. If the decay turns out to be exponential, it will be very intriguing because it will im- 
ply that an energy gap opens up owing to interactions. Note that a weak coupling analysis 
shows that the spectrum is continuous starting at zero even in a finite volume because of 
noncompact fiat directions of the classical potential energy.*^ Then this system provides an 
example in which the Witten index becomes well defined as a result of interactions while a 
perturbative analysis indicates that it is not so. 
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